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FIELD  STRENGTH  CALCULATIONS  FOR  E.  L.  F.    RADIO  WAVES 


Abstract 

The  mode  theory  of  propagation  of  electromagnetic  waves  at 
extremely  low  frequencies  (1.0  to  3000  c/s)is  briefly  reviewed  in  this 
paper.     Starting  with  the  representation  of  the  field  as  a  sum  of 
modes,   approximate  formulas  are  presented  for  the  attenuation  and 
phase  constants.     At  the  shorter  distances,    where  the  range  is 
comparable  to  the  wavelength,   the  spherical-earth  mode  series  is 
best  transformed  to  a  series  involving  cylindrical  wave  functions. 
This  latter  form  is  used  to  evaluate  the  near  field  behavior  of  the 
various  field  components.     The  effect  of  the  earth's  magnetic  field  is 
also  evaluated  using  a  quasi-longitudinal  approximation  which  is 
particularly  appropriate  for  propagation  in  the  magnetic  meridian. 
In  general  it  is  indicated  that  if  the  gyro  frequency  is  comparable  or 
greater  than  the  effective  value  of  the  collision  frequency,    the 
presence  of  the  earth's  magnetic  field  may  be  inaportant  for  E.L.F. 
In  this  case  the  attenuation  may  be  increased  somewhat.     The  in- 
fluence of  a  purely  transverse  magnetic  field  is  also  considered. 

This  technical  note  is  actually  a  numerical  supplement  to  the 
paper  "Mode  Theory  and  the  Propagation  of  ELF  Radio  Waves,  "  by 
J.R.    Wait,    J.    Research  N.  B .  S.  ,    64D,    ....    July/Aug.   I960;  however, 
for  sake  of  completeness  the  relevant  theory  is  briefly  presented. 

1.     Introduction 

The  propagation  of  electromagnetic  waves  at  extremely  low 
frequencies  (E.L.F.)  has  received  considerable  attention  recently 
[1-5]  .     Unfortunately,    the  experimental  data  in  this  frequency  range 
(1.0  c/s  to  3000  c/s)  is  still  rather  limited  [  8-15]  .     There  is  little 
doubt,   however,    that  the  bulk  of  the  energy  is  transferred  via  a 
waveguide  mode.     The  bounding  surfaces  of  this  waveguide  are  the 
ground  and  the  lower  edge  of  the  ionospheric  E  region. 
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An  important  feature  of  E.  L.  F.   propagation  is  that  the  distance 
between  source  and  observer  may  be  comparable  to  the  wavelength. 
For  example,   at  300  c/s  the  wavelength  is  1000  km.     In  fact,    many 
of  the  experiments  are  carried  out  at  distances  of  this  order  and 
account  must  be  taken  of  the  near  field  effects. 

It  is  the  purpose  of  this  note  to  present  computations  based  on 
the  mode  theory  as  it  applies  to  E.L.F.   propagation.     It  is  also 
hoped  that  these  results  will  also  be  of  assistance    in  interpretation 
of  future  experimental  data. 

II.    Basic  Theoretical  Model 

The  as  sunned  theoretical  model  is  taken  to  be  a  homogeneous 
conducting  spherical  earth  of  radius  a  surrounded  by  a  concentric 
conducting  ionospheric  shell  of  inner  radius  a  +  h.     It  is  convenient 
to  introduce  the  usual  spherical  coordinate  system  (r,  9 ,  (f>) .     Using 
this  model  and  assuming  a  vertical  electric  dipole  source  located  at 

Q  =  Q     r  =  a,   the  expression  for  the  vertical  electric  field,   E    ,    is 

r 

given  by  an  expression  of  the  following  form  [4,    16] 

E     =  14^    y       6      Mll±il    P     (-  COS  9)  (1) 

r       Zkha'^      /_,         n    sm  v  tr         v  '  ^  ' 

n  =  0 

where  I  -  -    average  current  in  the  source  dipole; 
ds  -  -  length  of  source  dipole; 

T|  -  -  intrinsic  impedance  of  air  space  =  IZOir  ohms; 
P  (-  cos  9)  -  -  Legendre  function  of  argument  -  cos  9  and 
(complex)  order  v; 


V  + -T-  =   ka  S_  where  S_,    for  n  =  0,   1,    2,    .    .    ,    , 


-  =   Ka  t  _    . 

n  n 

deternmined  from  the  boundary  condition  and  is 

described  below;  and 

^   1/2  for  n  =  0 

=  1  for  n   =    1,    2,    3  . 
'n 

The  individual  terms  in  the  above  series  correspond  to  the 

waveguide  modes.     In  the  general  case,   the  factor  S    is  obtained  from 


5 

n 

^               1 

sin  2  kh  C 
1  +                        ^ 

^       Z  kh  C 
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the  solution  of  a  complicated  transcendental  equation  which  involves 
spherical  Bessel  functions  of  large  argument  and  connplex  order. 
Certain  aspects  of  this  problem  have  been  discussed  by  Schumann  [  1] 
and  also  the  author  [  4]  .     For  a  homogeneous  earth  and  a  homogeneous 
ionosphere,   both  assumed  isotropic,    it  has  been  shown  that  S^  may  be 
approximated  by  [17] 


n 


s.-[-c^y^^^,K7l^^]y 


(2) 


where      A  =  ^+  ^ 
1  g 

in  terms  of  the  refractive  indices,    N    and  N.  ,    of  the  earth  and  the 

g  1 

ionosphere,    respectively.     This  equation  is  valid  subject  to  the  con- 
dition  j  a|   kh  «  1.      The  sign  of  the  radical  is  chosen  in  the  above 
equation  which  makes  the  real  part  positive.     The  values  of  S    then 
are  located  in  the  fourth  quadrant  of  the  complex  plane.     When  n  =  0, 
the  above  simplifies  to 


S     = 
o 


1 


.      A 


A 

"ZTOi 


(3) 


Now,    since    JA|  kh  «  1,    the  radical  in  equation  (1)  can  be  expanded 
for  n  >  0  to  yield 


S    = 
n 


1  - 


irn 
TEE 


2       .    2A 
-  1 


2 


er 


for  n  =  1,    2,    3 


(4) 


Furthermore,    if  in  addition, 


A 


«  1 


Trn  \ 
kE    ' 


it  is  permissible  to  write,    for  all  n, 

i  € 

2 


Sn" 


1  - 


Trn^' 
lEE 


n 


ZTcH 


/irn  \' 

(erj 


(5) 


where  €„  =  1.  e     =2  (for  n  =  1,    2,    3  .    .    .    ).     This  latter  equation  is 

valid  when  the  "waveguide  naodes"  are  not  near  cut-off  and  is  in  the 

form  given  by  Schumann [  2]  .     It  is  quite  analogous  to  the  standard 

results  for  the  propagation  constant  in  a  rectangular  waveguide  with 
finitely  conducting  walls  [18,    19]  . 
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III.     Earth-Flattening  Approximation 
The  possibility  that  the  curvature  of  the  earth  may  be  neglected 
is  now  investigated.     To  simplify  field  calculations  at  relatively  short 
distances,    it  would  seenn  desirable  to  transform  the  mode  series  to  a 
form  where  the  first  term  corresponds  to  the  model  of  a  flat  earth  and 
the  succeeding  ternas  are  corrections  for  curvature.     This  approach 
has  been  used  by  Pekeris  [20]    and  more  recently  by  Koo  and  Katzin 
[21]    in  their  investigations  of  microwave  duct  propagation.     From, 
their  work,   it  may  be  shown  that 

P    (-  cos  0)  1     /  o  \2      /    \ 

-^ =   H<2)^kS   p)   -    ^(-^)    H^2'(kSp)  (6) 

sm  V  IT  o     ^      n'^'         12  I    a    /        2  n^'  ^    ' 

4  ft 

+  terms  in  (p/a)   ,    (p/a)  ,  etc.  , 

where  v  +  l/l  =  kS  p  and  p    =   a0.     H^'  and  FT^'  are  Hankel  functions 
'  n^  ^  o  2 

of  the  second  kind  of  order  zero  and  two,     respectively.     When  the 
great  circle  distance  p  is  reasonably  small  compared  to  the  earth's 
radius  a,    only  the  first  term  in  the  expansion  need  be  retained. 

For  convenience  in  what  follows,    it  is    desirable    to  express  the 
field  components  as  a  ratio  to  the  quantity 


E 


o 


=    i  (nA)  I  ds  (e"'^P)/p,       (ti  =  12077); 


E    is  the  radiation  field  of  the  source  at  a  distance  p  on  a  perfectly 
conducting  ground.     Thus,    for  both  the  source  and  the  observer  near 
the  ground,    it  is  not  difficult  to  show  that 


E    =   W  E 
z  o  00 


where         W    =  -  itr^  e'^P  7   6    S^  H^2)^kS  p) 

h  i_j    n    n       o     ^     n^ 


(7) 


n=0 


E     =  -  S  E 
p  o 

OO 


where        S    ^   ^  ^  e^^P    ^    6^  S^  k}^^  (kS^p)  (8) 

^  n=0 

and  H .  =    -  T  E  /ti 

9  o'  ' 
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where 


T   ='    -  TT  f    e 


ikp 


00 


(2) 


>      5    S     H^2'    (kS  p)  =    -  N    S 


(9) 


n=o 


where  terms  containing  (p/a)"^  ,  (p/a)"*  have  been  neglected. 

When  kp  »  1,    corresponding  to  the  "far-zone,  "  the  above  ex- 
pressions may  be  simplified  since  the  Hankel  functions  may  be  re- 
placed by  the  first  term  of  their  asymptotic  expansion.     This  leads  to 
the  compact  result 

W 


S 
T 


^^'"^-^1      «n 


3/2 


n=o 


n 


n  '    g 
n      ___ 


-  iZirS  pA      (10) 


which  is  valid  for  p  »    \.     As  expected,   the  ratio  of  W  to  T  for  a 
given  mode  is  S    which  for  low  order  or  grazing  modes  is  of  the  order 
of  unity.     The  ratio  of  S  to  T,   quite  generally,    is   -  l/JST    which  is  very 
small  connpared  to  unity;  in  fact,    it  vanishes  for  a  perfectly  conduct- 
ing ground  as  it  must. 

The  excitation  by  a  horizontal  electric  dipole  may  be  treated  in 
a  similar  manner.     A  complete  discussion  of  the  theory  is  given  else- 
where [17]  .     The  expression  for  the  vertical  field  E    of  the  horizontal 
dipole  Ids  which  is  located  at  the  origin  and  oriented  in  direction 
0  =  0  may  be  written  in  the  form 

E      =    S  E    cos  0  (11) 

z  o 

where  S  and  E    are  defined  above.     By  using  the  reciprocity  theorem, 
this  result  can  also  be  deduced  directly  from  equation  (9).     Compar- 
ing the  above  formula  for  E    (for  a  horizontal  dipole)  with  correspond- 

z 

ing  formula  for  E  ,    for  a  vertical  dipole,    it  nnay  be  easily  shown  that 


E    (for  horizontal  dipole) 
z 

E    (for  vertical  dipole) 


cos 


N 


0     /  1   \2  cos 

~  ( ^  "1^/   ^  "n: 


S    0 


(12) 


z  g        ^       '\  '  g 

being  valid  in  the  asymptotic  or  far  zone.     In  terms  of  the  ground 
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c onductivity  <r     and  dielectric  constant  c    , 

o  o 

g  g       g  -     g 

This  is  a  small  quantity. 

The  above  formulas  for  W,   S,   and  T  are  valid  only  if  p  «  a. 
If  the  first  curvature  correction  term  is  included,    it  is  a  simple 
matter  to  show,    in  the  "far  zone",   that  this  amounts  to  naultiplying 
the  right-hand  side  by  the  factor 

■  ^  / 

\ 

Starting  from  the  mode  series  for  a  spherical  earth  [  17]    and  using 
the  far  field  approximation  for  the  Legendre  function,   it  turns  out 
that  W,   S,    T  have  precisely  the  sanne  fornn  as  in  equation  (10)  ex- 
cept the  factor  {9/sin9)^  occurs  in  the  right-hand  side.     Noting  that 
1 


2 

terms  in  ^,    0"  ,    etc, 


ey  .!+«!+•         •    -^  -^ 


"sHT^  ~      ^  TT 

it  is  apparent  to  a  second  order  that  the  derived  curvature  correction 
is  the  same  as  from  the  Hankel  function  series. 

IV.     Distance  and  Frequency  Dependence 

Certain  features  of  the  mode  series  are  best  illustrated  by  a 
graphical  plot  of  field  strength  vs.    distance  under  various  values  of 
the  parameters.     Since  N  »  N.  it  may  be  safely  assumed  for  E.  L.  F. 
that  A  =  1/N.   •       Then  on  the  assumption  that  the  ionosphere  is  be- 
having like  a  conductor,   it  is  possible  to  write 


-i  -^^ 


where 

_  cOq^  _   (plasma  frequency)^ 
r  V  collision  frequency 
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The  above  relation  for  N.   is  an  approximation  which  is  valid  when  the 

collision  frequency  v  between  electrons  and  neutral  ions  is  much  less 

7 
than  the  angular  frequency  o).     Since  v  is  of  the  order  of  10     in  the 

lower  E  region,    this  is  certainly  valid  at  E.L.F.      The  (angular) 

plasma  frequency  co,,   is  given  by 

2  e^N 


00, 


€o   m 


where  N  is  the  electron  density,    e  and  m  are  charge  and  m.ass  of  the 
electron,    and€„   is  the  dielectric  constant  of  free  space,     oj   ,    which 
is  proportional  to  the  ratio  N/v,    is  believed  to  be  of  the  order  of  10 
in  the  D  layer  or  the  lower  region  of  the  E  layer.     At  least  this 
appears  to  be  the  effective  value  deduced  from  V.  L.  F.    observations 
for  highly  oblique  incidence  [  Z3]  .     At  E.  L.F.    it  appears  that  the 
waves  are  reflected  at  higher  levels  in  the  ionosphere  and  the  effect- 
ive value  of  CO    is  sonaewhat  larger,   being  possibly  of  the  order  of 
10^ 

It  should  be  noted  in  passing  that  co    may  be  replaced  by  (r-Ao 
where  q-.  is  the  effective  conductivity  of  the  ionosphere  and  ^   is  the 
dielectric  constant  of  free  space. 

To  indicate  the  variation  of  the  field  as  a  function  of  distance, 
the  magnitude  and  phase  of  T  and  W  are  shown  in  Figs,    la  to  4c  for  a 
range  of  values  of  co    and  h  =  90  km.     The  curves  denoted  co    =  oo 
correspond  to  a  perfectly  conducting  ionosphere.      These  are 

presented  for  purposes  of  comparison  with  the  set  for  the  more 

6  5 

reasonable  values  of  co  .     These  are  2  x  10     and  5  x  10    .     Actually, 

|t|   and   |wj   are  both  divided  by  the  square  root  of  the  distance  p 

(in  kilometers)  and  plotted  since  the  curves  become  linear  at 

larger  distances.      The  slope  of  these  linear  parts  of  the  curves  is 

proportional  to  the  attenuation  rate  (in  db  per  1000  km)  as  usually 

defined.     It  is  seen  immediately  that  at  the  shorter  distances,    the 

curves  are  no  longer  straight.     This  immediately  indicates  that 

considerable  caution  should  be  exercised  in  computing  attenuation 

rates  from  spectral  analyses  of  "sferics.  "    For  example,    it  is  only 
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when  the  distance  p  exceeds  about  one -sixth  of  a  wavelength,    is  it 

permissible  to  as surae  that  log     (|w|/'"Tp~)    or  log     (Ie   |  x  's/p  )  vary 

z 

in  a  linear  manner  with  distance  p.     Similarly,   the  phase  of  W(orE  ) 
varies  in  a  linear  manner  only  when  the  distance  exceeds  about  one- 
half  wavelength.     Similar  remarks  apply  to  the  magnetic  field 
variations  with  distance. 

To  shed  further  light  on  the  nature  of  the  E.L.F.    fields  in  the 
waveguide,   the  radial  impedance  of  the  wave  is  now  considered. 
It  is  noted  that 


E 


=   -Z  (14) 

z  =  0  ^ 


which  is  by  definition  the  (normalized)  impedance  of  the  wave  looking 
in  the  radial  or  p  direction.     First  it  should  be  noted  that  for  kp  »   1, 

^^  1 

T 
and  for  kp  «   1  and  p  «  h, 

W    ^      1 


T  ikp 

For  the  intermediate  and  interesting  range,   kp  is  comparable  to  unity 
and  p  is  somewhat  greater  than  h.     Using  the  numerical  values  of  T 
and  W  nnentioned  above,   the  ratio   |  w/ T  |    and  the  phase  (lag)  defined 
by  arg  T  -  arg  W  are  plotted  in  Figs.    5a  to  6c  and  as  a  function  of 
distance  for  various  frequencies.     As  before,   h  =  90  km  and  the  same 
range  of  co    values  are  chosen. 

The  impedance  ratio  W/T  is  independent  of  the  frequency 
spectrum  of  the  source  provided,    of  course,   it  may  be  represented 
by  an  equivalent  vertical  electric  dipole.     In  fact,   this  complex  ratio 
could  be  easily  calculated  from  the  frequency  spectra  of  the  wavefornas 
of  the  vertical  electric  field  and  the  horizontal  magnetic  field  of  an 
atmospheric.     The  observed  variation  of  magnitude  or  phase  of  W/T 
as  a  function  of  frequency  should  then  provide  a  basis  for  distance 
measuring.     Such  a  scheme,   while  admittedly  crude,    requires  only 
one  receiving  station  equipped  with  a  vertical  whip  and  a  loop  antenna. 
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At  the  very  short  distances  where  p  is  of  the  order  of  the 

ionospheric  reflecting  height  (i.e.  ,    90  km),    the  mode  series 

representation  for  the  fields  is  very  poorly  convergent.     However, 

alternate  representations  are  available  which  converge  very  rapidly 

in  this  case  [  17]  .     For  co    =   oo,   these  are 


T    = 


m 


m=o 


_t^Vl-^^e-i'^<'■m-P' 


m 


(15) 


m 


and 


00 


W    = 


m=o 


L  r       ) 

[(■  - .: 

1 

\        i 

■ ; "  RF- 

m  J 

/2mhV 
\      m 

[6- 

3i 

"    kr 
m 

'p^+( 

2mh)2 

1 

=  1,     e        = 
m 

=  2(m 

(16) 


ik(r      -p) 


m 


where  r      = 
m 


Individual  terms  in  these  series  correspond  to  the  contribution  from 

images  of  the  source  dipole.     These  particular  expressions  are  exact 

on  the  assumption  of  perfectly  conducting  planes,      (i.e.  ,    oj    =   oo  and 

cr     =  oo).    It  is  known,   however,    from  the  mode  calculations  that  at 

g 
short  distances  (i.e.  ,    for  p  <   JOO  km),    the  numerical  values  of  the 

functions  T  and  W  are  only  very  slightly  dependent  on  the  finite  value 

of  u)    and  o"    .       Thus  (15)  and  (16)  are  generally  applicable. 

To  illustrate  the  nature  of  the  fields  at  short  ranges,    the  magni- 
tude and  phase  of  T  and  W  are  shown  plotted  in  Figs.    7a  to  8b.     The 
solid  curves  correspond  to  equations  (15)  and  (16)  for  h  =  90  km.     The 
dashed  curves  correspond  to  equations  (15)  and  (16)  when  h  =  oo;  this 
is  the  sanae  as  taking  only  the  first  term  in  the  series.     At  the  very 
short  ranges,    it  is  seen  that  the  solid  and  dashed  curves  merge  to- 
gether.    Apparently  in  this  region,    the  reflecting  layer  at  90  km  can 
be  ignored.     At  the  greater  ranges,    the  T  and  W  functions  both 
approach  unity  for  the  case  of  no  reflecting  layer  (i.e.,   h  =  oo). 
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Corrections  for  earth  curvature  resulting  from  diffraction  would, 
of  course,    modify  these  dashed  curves  at  ranges  exceeding  about 
500  km.    It  may  be  concluded  from  the  information  plotted  in  Figs. 
7  and  8  that  only  for  extremely  short  ranges  (i.e.  ,    less  than  about 
30  km)  is  it  permissible  to  neglect  the  presence  of  the  ionospheric 
reflecting  layer  in  field  strength  calculations  at  ELF. 

V.     Effect  of  the  Earth's  Magnetic  Field 

In  the  preceding  it  has  been  tacitly  assumed  that  the  ionosphere 
is  behaving  as  an  isotropic  homogeneous  conductor.     In  this  case, 
the  refractive  index  N.  of  the  ionosphere  may  be  written  [23] 

1 


N.     s 


1  i. 

i  _£.  I        s  /  .JL  '  (17) 

0)  I      10) 


where 

_    cof      _    (Plasma  Frequency)^ 
r  V  (Collision  Frequency) 

When  the  earth's  magnetic  field  is  steeply  dipping,    it  is  appropriate 

to  invoke  the  quasi-longitudinal  approximation  of  Booker  [23]  .     In 

this  case  the  refractive  index  is  double  valued,    one  corresponds  to 

the  ordinary  and  the  other  the  extra-ordinary,    thus 


Nj^ 


r 
1  -  i    exp  (±  It) 


O) 


(18) 


where 

tan  T  = 

and 


L,       _     longitudinal  connponent  of  gyro -frequency 
V        ~  collision  frequency 


CM? 


O       = j_       . 

^  (v2+cu|^)2 

Using  this  model  for  the  ionosphere,  the  reflection  coefficient  for  a 
sharp  boundary  has  been  derived  by  Budden  [24]  .  Adapting  this  re- 
sult to  E.Li.F.   it  has  been  shown  [17]    that  the  formula  for  S    has  the 

^     ^  n 
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same  form  as  the  isotropic  medium  if  A  is  defined  by 

_i 

g  V  "r 

Since  O      =    u)     cos  t  it  is  possible  to  write  the  results  in  the  same 
r  r  ^ 

form  as  the  isotropic  ionosphere  if  we  set 


where 


(Ni)e£f.     = 

^    r'  efi.      ^ 
ico 

(21) 


(N.)    rr    and  (u)  )    ,^    are  the  effective  values  of  N.  and  co   ,    respectively. 
*    reff.  r'eff.  i  r  ^  •' 

Specifically, 

^    r'eff-  r     (cos  t/2)2  ' 

In  other  words,   the  effective  conductivity  of  the  ionosphere  is  modi- 

COS    T 

fied  by  the  factor  -, /-,  ^„     which  varies  fronn.  unity  to  zero  as  t 

^  (cos  t/2  )2  -' 

varies  from  0  to   ir/Z.     Thus, the  attenuation    is  increased  as  a  result 
of  a  steeply  dipping  (or  vertical)  magnetic  field.     In  fact 

Attenuation  with  magnetic  field  ^      cos  (t/Z) 

Attenuation  without  nnagnetic  field      ~      ,  ,1 

^  (cos    T  )   2 

If  T  «   1  (i.  e.  ,    (ii     «   v),   this  ratio  becomes  unity  and  the  influence 

of  the  earth's  magnetic  field  vanishes.     At  the  level  in  the  ionosphere 

where  E.  L.  F.    waves  are  reflected,    it  is  expected  that  co     and  v  are 

comparable,   both  being  of  the  order  of  10    .     Assuming  that  they  were 

actually  equal,    t  becomes  4  5     and  the  ratio  of  the  attenuation  rates 

is  1.05.     Thus,    it  is  only  when  o)      is  greater  than  v  does  the  earth's 

J— I 

magnetic  field  appreciably  influence  the  attenuation.     The  ratio  is 
plotted  in  Fig.  9    as  a  function  of  tan  t  or  oj    /v.     This  illustrates  the 
situation  clearly. 

The  preceding  results  are  subject  to  the  validity  of  Booker's 
quasi -longitudinal  approximation  of  the  Appleton-Hartree  equation[23] 
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The  validity  of  this  approximation  requires  that 
4 


where  cj^    and  co     are  the  longitudinal  and  transverse  components  of 
the  (angular)  gyro-frequency.     Clearly  this  condition  is  violated  when 
the  transverse  component  of  the  earth's  nnagnetic  field  is  large  such 
as  for  propagation  along  the  magnetic  equator.     This  case,   however, 
has  been  considerrd  by  Barber  and  Crombie[25]    who  derived  explicit 
results  for  the  reflection  coefficient  at  a  sharply  bounded  ionosphere 
with  a  purely  transverse  magnetic  field.     Adapting  their  results  to 
E.  L.  F.    it  is  not  difficult  to  show  [17]    that  A  is  now  given  by 

A=g)*     X  (23) 

CU/l  CO  v2co/  VVW 


X    ^ 


J       \  r^\  r  r^  \r 


2 


\   2  '^rp  2 

,    ,    .     CO     \  T     oj 

1  +  1  —       _ 

r    '  r 

The  real  and  imaginary  parts  of  A  are  plotted  in  Figs.  10a  and  10b  as 
a  function  of  the  frequency  parameter  5x  10    x  f/co     where  f  is  ex- 
pressed in  kc/s  and  co     is  in  (sees)      .     The  various  values  of  a)_/v  are 

r  X 

shown  in  the  curves;  positive  values  of  the  ratio  correspond  to 
propagation  from  east-to-west  while  negative  values  correspond  to 
propagation  from  west-to-east.     Now  the  attenuation  of  the  dominant 
mode  in  nepers  per  unit  distances  is  given  by  -  kS    where  S    is  given 
by  equation  (3). 

Using  the  above  results  for  A  the  attenuation  coefficient  q,    ex- 
pressed in  db  per  1000  km  of  path  length,    is  plotted  in  Fig.    Ha  for 
h  =  90  km.     It  is  immediately  evident  that  the  attenuation  is  only 
slightly  raodified  by  the  presence  of  a  transverse  magnetic  field 

unless,    of  course,  co_p  is  considerably  greater  than  v.     The  phase 

T 

velocity  of  the  dominant  zero-order  mode,    relative  to  free  space,    is 
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simply  l/Re    S    .     Denoting  this  dinriensionless  quantity  by  p  it  is 
plotted  in  Fig. lib    for  h=  90  km.     It  is  apparent  that  the  presence  of 
a  transverse  magnetic  field  does  have  an  influence  on  the  phase 
velocity. 
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Figure  Captions 

Figs,   la,  b,   c.        The  normalized  magnitude  of  the  magnetic  field  as 

a  function  of  distance  from  the  source  (p  is  ex- 
pressed in  kilometers  in  both  ordinate  and  abscissa) 

Figs.    2a,  b,    c,       The  normalized  phase  (lag),  in  radians,    of  the 

magnetic  field. 

Figs.    3a,   b,    c,       The  normalized  magnitude  of  the  electric  field  as  a 

function  of  distance  from  the  source. 

Figs.   4a,   b,    c,       The  normalized  phase  (lag)  of  the  electric  field. 
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Figs.    5a,   b,    c         The  magnitude  of  the  radial  wave  impedance. 

Figs.    6a,   b,    c         The  phase  (lag)  of  the   radial  wave  impedance. 

Figs.    7a,   b  The  magnitude  and  phase  of  the  normialized  magnetic 

field  at  short  distances.     The  broken  curves  corres- 
pond to  no  reflecting  layer. 

Figs.   8a,   b  The  magnitude  and  phase  of  the  normalized  electric 

field  at  short  distances. 
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Fig.    "^  The  earth's  -magnetic -field-correction  to  the  attenuation 

factor  for  propagation  in  the  magnetic  meridian. 
Figs.    10a.   b    The  real  and  imaginery  parts  of  the  factor  A  as  a 

function  of  normalized  frequency  (for  a  purely  transverse 

magnetic  field). 
Figs.  11a,   b      Attenuation  and  relative  phase  velocity  for  a  purely 

transverse  field. 
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Atomic  and  Radiation  Physics.  Spectroscopy.  Radiometry.  Mass  Spectrometry.  Solid  State  Physics. 
Electron  Physics.  Atomic  Physics.  Neutron  Physics.  Radiation  Theory.  Radioactivity.  X-rays,  fligli 
I  Energy  Radiation.     Nucleonic  Instrumentation.     Radiological  Equipment. 

Chemistry.  Organic  Coatings.  Surface  Chemistry.  Organic  Chemistry.  Analytical  Chemistry.  Inorganic 
Chemistry.  Electrodeposition.  Molecular  Structure  and  Properties  of  Gases.  Physical  Chemistry.  Thermo- 
chemistry.    Spectrochemistry.     Pure  Substances. 

Mechanics.  Sound.  Mechanical  Instruments.  Fluid  Mechanics.  Engineering  Mechanics.  Mass  and  Scale. 
Capacity,  Density,  and  Fluid  Meters.     Combustion  Controls. 

Organic  and  Fibrous  Materials.  Rubber.  TextUes.  Paper.  I,eather.  Testing  and  Specifications.  Poly- 
mer Structure.     Plastics.     Dental  Research. 

-Metallurgy.  Thermal  Metallurgy.  Chemical  Metallurgy.  Mechanical  Metallurgy.  Corrosion.  Metal 
Physics. 

Mineral  Products.  Engineering  Ceramics.  Glass.  Refractories.  Enameled  Metals.  Constitution  and  Mi- 
crostructure. 

Building  Technology.  Structural  Engineering.  Fire  Protection.  Air  Conditioning,  Heating,  and  Refrigera- 
tion.    Floor,  Roof,  and  Wall  Coverings.     Codes  and  Safety  Standards.     Heat  Transfer.     Concreting  Materials. 

Applied  Mathematics.     Numerical  Analysis.     Computation.     Statistical  Engineering.     Mathematical  Physics. 

Data  Processing  Systems.  SEAC  Engineering  Group.  Components  and  Techniques.  Digital  Circuitry. 
Digital  Systems.     Analog  Systems.     Application  Engineering. 

.  Office  of  Basic  Instrumentation.  ,  office  of  Weights  and  Measures. 
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Cryogenic  Engineering.  Cryogenic  Equipment.  Cryogenic  Processes.  Properties  of  Materials.  Gas  Lique- 
faction. 

Radio  Propagation  Physics.  Upj)er  Atmosphere  Research.  Ionospheric  Research.  Regular  Propagation 
Services.  Sun-Earth  Relationships.  VHF  Research.  Radio  Warning  Services.  Airglow  and  Aurora.  Radio 
Astronomy  and  Arctic  Propagation. 

Radio  Propagation  Engineering.  Data  Reduction  Instrumentation.  Modulation  Research.  Radio  Noise. 
Tropospheric  Measurements.  Tropospheric  Analysis.  Propagation  Obstacles  Engineering.  Radio-Meteor- 
ology.    Lower  Atmosphere  Physics. 

Radio  Standards.  High  Frequency  Electrical  Standards.  Radio  Broadcast  Service.  High  Frequency  Im- 
pedance Standards.     Electronic  Calibration  Center.     Microwave  Physics.     Microwave  Circuit  Standards. 

Radio  Communication  and  Systems.  Low  Frequency  and  Very  Low  Frequency  Research.  High  Fre- 
quency and  Very  High  Frequency  Research.  Ultra  High  Frequency  and  Super  H  igh  Frequency  Research.  Modula- 
tion Research.     Antenna  Research.     Navigation  Systems.     Systems  Analysis.     Field  Operations. 


